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Detonation  initiation  on  the 
microsecond  time  scale:  DDTs 

D.  R.  KASSOY*,  J.  A.  KUEHN*,  M.W.  NABITY',  J.  F.  CLARKE* 

Spatially  resolved,  thermal  power  deposition  of  limited  duration  into  a  finite  volume 
of  reactive  gas  is  the  initiator  for  a  deflagration-to-detonation  transition  (DDT)  on  the 
microsecond  time  scale.  The  reactive  Euler  equations  with  one-step  Arrhenius  kinetics 
are  used  to  derive  novel  formulas  for  velocity  and  temperature  variation  that  describe  the 
physical  phenomena  characteristic  of  DDTs.  A  transformation  of  the  variables  is  shown 
to  yield  a  canonical  equation  system,  independent  of  the  activation  energy.  Numerical 
solutions  of  the  reactive  Euler  equations  are  used  to  describe  the  detailed  sequence  of 
reactive  gasdynamic  processes  leading  to  an  overdriven  planar  detonation  far  from  the 
power  deposition  location.  Results  are  presented  for  deposition  into  a  region  isolated 
from  the  planar  boundary  of  the  reactive  gas  as  well  as  for  that  adjacent  to  the  boundary. 
The  role  of  compressions  and  shocks  reflected  from  the  boundary  into  the  partially 
reacted  hot  gas  is  described.  The  quantitative  dependences  of  DDT  evolution  on  the 
magnitude  of  thermal  power  deposition  and  activation  energy  arc  identified. 
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1  Introduction 

“Gasdynamics  of  explosions  is... best  defined  as  the  science  dealing  with  the 
interrelationship  between  energy  transfer  occurring  at  a  high  rate  in  a  compressible 
medium  and  the  concomitant  motion  set  up  in  the  this  medium"  (Oppenheim  and 
Soloukhin.  1973).  The  dynamics  of  the  interaction  depend  on  the  relationship  between 
the  dimensional  time  scale  for  substantial  chemical  heat  release,  t'H ,  into  a  local  region  of 
length  scale  f ,  and  local  acoustic  time  l'A  =Cja  where  a  is  the  local  speed  of  sound. 
When  t'H  <zt'A  local  heat  addition  occurs  in  a  basically  constant  volume  process  (local 

inertial  confinement)  because  the  fluid  cannot  respond  on  the  short  time  scale  (Kassoy 
and  Palaniswamy,  2002).  The  temperature  rise  is  accompanied  by  a  concomitant 
pressure  rise,  so  that  for  a  brief  instant  a  high-pressure  spot  exists  in  a  relatively  low 
pressure  and  temperature  environment.  Subsequent  expansion  of  the  spot  on  the  local 
r'^-time  scale,  driven  by  the  large  pressure  gradient  between  the  spot  and  the  environment 

is  the  source  (“piston"  effect)  of  compression  waves  in  the  environment.  Wave 
coalescence  can  lead  to  shock  wave  formation.  Partial  inertial  confinement  occurs  in  the 
less  exteme  case  l'tl  =  0(iA )  so  that  localized  pressure  increase  with  heat  release  is  still 

possible  although  the  expansion  process  occurs  simultaneously  with  the  heat  addition, 
limiting  the  pressure  rise  and  the  amplitude  of  the  compression  waves.  “Concomitant  to 
(these  ideas)  is  the  appreciation  of  the  role  played  in  the  detonation  process  by  the  power 
density.. ..at  which  energy  is  deposited  into  the  compressible  medium”  (Oppenheim  and 
Soloukhin,  1973).  This  quote  refers  to  reactive  gasdynamic  transients  specifically, 
where  gasdynamic  w'ave  generation  following  a  high  rate  of  localized  heal  addition  is 
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crucial  to  facilitating  fast,  high  temperature  reaction  rates,  in  contrast,  the  chemical  heat 
release  rate  is  insignificant  for  a  steadily  propagating  Chapman- Jouguct  (C.J.)  wave 
((■‘ickelt  and  Davis,  1979). 

The  evolutionary  process  of  detonation  formation  following  thermal  power 
deposition  into  a  reactive  mixture  is  fundamentally  dependent  on  reactive  gasdynamic 
transients.  Experiments  of  this  type  in  a  tube  (Oppenheim.1985)  begin  with  a  laminar 
flame  initiated  by  a  spark.  A  myriad  of  instabilities  cause  the  flame  to  become  wrinkled 
(Khokhlov,  Oran,  and  Thomas,  1999,  and  Khokhlov  and  Oran.  1999),  thus  enhancing  the 
rate  of  energy  production  through  area  enhancement.  The  resulting  turbulent  flame  brush 
is  observed  to  be  the  source  of  gasdynamic  compression  waves  propagating  forward  into 
the  initially  undisturbed  reactive  mixture.  The  local  reaction  rate  is  enhanced  by 
temperature  increases  caused  by  compression  wave  passage.  Eventually  a  definitive 
localized  “explosion  in  an  explosion"  occurring  on  the  tube  boundary  is  the  source  of 
strong  shocks  which  coalesce  with  previously  generated  compression  waves  to  form  a 
shock  capable  of  inducing  substantial  chemical  heat  release  immediately  downstream  of 
the  wave  front.  It  is  possible  that  a  transient  ZND  detonation  wave  is  produced  first. 
However,  in  experiments  all  such  waves  evolve  to  a  three  dimensional  cellular  structure, 
Lee,(l984).  A  very  considerable  literature  describes  experimental  and  modeling  studies 
of  one  and  two  dimensional  instabilities  of  initially  planar  ZND-waves.(see,  for  example, 
Gamezo,  Desbordes  and  Oran,  1999.  Sharpe  and  Falle.  1999,  2001,  Sharpe,2001, 
Kasimov  and  Stewart,  2004  for  reviews  of  major  contributions).  The  emphasis  is  on 
highly  accurate  algorithms,  capable  of  resolving  the  evoution  of  a  steady  C.J.  wave  to 
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either  one  dimensional  transient  dynamics  predicted  by  linear  stability  theory(  Lee  and 
Stewart.  1990,  Sharpe,  1997,  Stewart  and  Kasimov,  2006)  and  to  two  dimensional 
transient  dynamics,  including  triple  points  and  other  transverse  wave  phenomena. 
Stability-related  modeling  is  distinct  from  that  describing  the  appearance  of  planar  or  two 
dimensional  detonations  following  a  defined  localized  thermal  disturbance,  of  interest  for 
facilitating  robust,  reliable,  repetitive  detonation  initiations  in  a  pulsed  detonation  engine. 
(Schauer.  Miser,  Tucker,  Bradley  and  Hoke.  2005). 

Dold  et  al.  (1995)  have  used  an  initial-boundary  value  solution  approach  for  the 
Euler  equations  to  find  one-dimensional  detonations  arising  from  an  initially  imposed, 
spatially  distributed  pressure  pulse.  Although  thermal  initiation  is  absent,  the  imposed 
pressure  disturbance  steepens  into  shocks  that  create  the  high  temperature  regimes 
necessary  for  rapid  acceleration  of  reaction  rates  and  the  appearance  of  multiple  ignition 
sites.  The  sequence  of  events  leading  to  detonation  formation  are  more  than  reminiscent 
of  the  events  in  thermal  initiation  models  by  Clarke,  Kassoy,  and  Riley  (1986)  and 
Clarke,  Kassoy,  Meharzi,  Riley,  and  Vasantha  (1990),  and  Sileem.  Kassoy,  and  Hayashi 
(1991)  (SKH).  This  qualitative  agreement,  with  in  fact  quite  similar  time  scales  of 
transition  to  detonation  for  similar  activation  energy  values,  is  very  encouraging  because 
the  numerical  method  employed  is  unlike  the  MacCormack  scheme  used  in  SKH. 

Sileem,  Kassoy  and  Hayashi  (I991)(SKH)  model  the  evolution  of  planar  detonation 
formation  following  time  and  spatially  resolved  thermal  power  into  a  finite  volume  of 
reactive  gas.  The  reactive  Euler  equations  with  one-step  Arrhenius  kinetics,  scaled  so 
that  t\  =0(t\)  are  solved  numerically.  The  evolution  to  detonation  is  described  in  terms 


4 


of  a  sequence  of  localized  chemical  heat  addition  events  leading  to  gasdynamic  wave 
generation  and  strong  shock  formation,  reminiscent  of  Oppenheim’s(  1985)  observations. 
An  overdriven  ZND-detonation  is  born  initially,  similar  to  the  experimental  observations 
of  Moen,  Donato,  Knystautas  and  Lee  (1980)  and  Schauer  et.  al.  (2005),  in  the  context  of 
an  evolving  turbulent  flame.  The  transient  ZND  wave  relaxes  spontaneously  to  a  stable 
CJ  wave  appropriate  to  the  mixture  parameters.  These  results  demonstrate  that,  from  a 
mathematical  perspective,  transverse  waves  and  turbulence  arc  not  essential  for 
detonation  formation.  A  physically  reasonable  sequence  of  planar  reactive  gasdynamic 
processes  featuring  accelerating  localized  chemical  power  deposition,  concomitant 
pressure  increases  and  subsequent  compression  wave  generation,  is  sufficient  to  ensure 
the  outcome. 

Smirnov  and  Panfilov  (1995)  describe  planar  deflagration  to  detonation  transition  in 
a  chemically  reactive,  viscous,  conductive  gas  with  a  two  step  kinetic  model.  Thermal 
initiators  included  boundary  heating,  like  that  in  Clarke  et  al.  (1986,  1990)  and  a  pre¬ 
existing  localized  high  temperature  spot.  A  pre-existing,  small,  localized  temperature 
gradient  is  used  by  Kapila  et  al.  (2002)  to  initiate  a  planar  detonation,  as  first  suggested 
by  Zeldovich  (1970,  1980).  Euler  equations  and  a  one-step  kinetic  model  are  used.  The 
Kapila  et  al.  paper  includes  an  extensive  literature  survey  of  related  work  for  both 
deflagration-to  detonation  (DDT)  processes  and  direct  initiation.  More  recently,  Gu  et  al. 
(2003)  have  modeled  the  formation  of  spherical  detonations  arising  from  a  pre-existing 
localized  hot  spot  on  the  millimeter  length  scale.  The  initial  spot  temperature  distribution 
is  characterized  by  a  linear  temperature  gradient  with  a  maximum  temperature  only  very 
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slightly  higher  than  that  of  the  surrounding  environment.  The  chemically  reactive, 
viscous,  conductive  gas  equations  are  used  along  with  a  hydrogen-oxygen  kinetic 
scheme.  Chemical  heat  release  occurs  on  a  microsecond  time  scale.  The  evolutionary 
process  to  detonation  is  found  to  be  quite  sensitive  to  the  characteristic  initial  gradient 
present. 

The  two-dimensional  reactive  Navier-Stokes  equations  with  one  step  kinetics  are 
used  by  Tegner  and  Sjogreen  to  model  transitions  from  deflagration  to  detonation  arising 
from  a  suddenly  imposed  hot  spot.  An  “accelerating  (laminar)  flame'’  is  the  initial  source 
of  a  lead  shock  wave  which  “preconditions”  the  unburned  mixture.  The  aforementioned 
reactive  gasdynamics  processes  follow. 

Sharpe  and  Falle  (1999)  describe  one-dimensional  evolution  of  an  initially  CJ  wave, 
based  on  solutions  to  the  Euler  equations  with  one  step  chemistry.  Numerical  noise 
provides  the  disturbances  that  yield  planar  transient  responses  reminiscent  of  those  seen 
in  the  SKI  1  initiation  study,  including  the  appearance  of  unburned  pockets  of  reactive 
mixture  well  behind  the  gasdynamic  wave  front.  The  irregularity  of  the  transients, 
defined  by  the  time  history  of  the  post  shock  pressure,  is  strongly  dependent  on  the 
activation  energy.  These  one-dimensional  simulations  require  significant  reaction  zone 
resolution. 

Kasimov  and  Stewart  (2004)  use  a  novel  numerical  algorithm  for  integration  of  the 
reactive  one-dimensional  Euler  equations  with  one  step  kinetics  in  a  shock-attached 
frame,  to  describe  “....the  dynamics  of  self-sustained  detonation  waves...."  Conditions 
on  the  upstream  boundary  are  given  careful  attention  to  minimize  the  impact  of  spurious 
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reflected  gasdyatnic  waves  on  the  interior  solution.  Solution  accuracy  requires  several 
hundred  grid  points  per  half-reaction  length.  Higher  activation  energy  mixtures  exhibit 
irregular  reactive  gasdynamics.  So  called  “re-ignition"  events  are  characterized  by  the 
appearance  of  localized  high  pressure  spots,  shown  in  big.  9,  which  are  correlated  with 
locally  high  rates  of  chemical  heat  release. 

Pre-existing  strong  blast  waves  can  be  employed  to  create  detonations  in  a  more 
immediate  way  than  is  typical  in  a  DDT.  Mazaheri  (1997)  and  Eckett  et  al.  (2000)  have 
modeled  direct  initiation  of  planar  and  spherical  detonations,  respectively,  following  the 
instantaneous  appearance  of  a  blast  wave.  I  he  modeling  is  based  on  one-step  kinetics. 
These  studies  predict  critical  energies  required  for  seif-sustained  detonations  as  well  as 
evolution  scenarios  for  failing  detonations, 

There  is  an  extensive  related  literature  describing  initial-boundary  value  solutions  to 
the  Euler  equations  for  mechanically  driven  plane  detonations.  Typically,  a  piston 
motion  is  prescribed  in  order  to  propagate  a  relatively  strong  shock  into  the  reactive  gas. 
Singh  and  Clarke  (1991)  provide  a  review  of  the  subject,  as  well  as  numerical  solutions 
found  by  using  a  random  choice  method. 

Related  efforts  in  computational  modeling  of  two-dimensional  initiation  of,  and 
transition  to  detonation,  are  described  by  Nikiforakis  and  Clarke  (1996).  A  computational 
solution  is  obtained  to  an  initial-boundary  value  problem  described  by  the  reactive  Euler 
equations.  The  initiator  in  this  case  is  a  planar  reflected  shock  propagating  through  a  spot 
with  an  imposed,  slightly  elevated  temperature  located  adjacent  to  the  closed  endwall  of  a 
shock  tube  in  an  asymmetric  position,  The  relatively  high,  spatially  non-uniform 
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temperature  behind  the  reflected  shock  accelerates  the  local  reactions  and  thereby 
initiates  the  desired  transient  reactive  gasdynamic  events.  These  planar  studies  predict 
wave  patterns  like  those  seen  in  the  observations  of  Meyer  and  Oppenheim  (1971 ). 

Gamezo  et.al.  (1999)  use  the  Euler  equations  with  one  step  chemistry  to  model  the 
evolution  of  an  initially  overdriven  planar  detonation  subjected  to  numerical  noise.  After 
a  period  of  attenuation  toward  a  steady,  planar  CJ  wave,  two-dimensional  perturbations 
of  the  shock  front  and  reaction  zone  begin  to  appear.  The  detonation  front  evolves  to  an 
unsteady  ensemble  of  transverse  waves  and  triple-shock  configurations  propagating  at  an 
average  axial  speed  of  a  planar  CJ  speed.  The  regularity  of  the  cell  structure  and  the 
transient  characteristics  of  the  transverse  phenomena  on  the  wave  front  are  quite  sensitive 
to  the  activation  energy. 

Strong  transverse  wave  evolution  from  a  disturbed  planar  CJ  wave  is  also  considered 
by  Sharpe  (2001).  The  two  dimensional  reactive  Euler  equations  with  a  one  step  kinetics 
model  are  solved  numerically.  A  spatially  resolved  density  inhomogeneity  is  used  to 
disturb  the  initially  planar  wave.  High  resolution  simulations  are  used  to  examine  the 
two  dimensional  structure  of  the  wave  front,  including  collisions  of  the  triple  points 
propagating  transversely,  and  transverse  phenomena  just  behind  the  front. 

The  primary  objective  of  the  current  modeling  study  is  to  extend  the  SKH  study  to  a 
broader  range  of  parameters,  including  the  initiating  thermal  power  deposition,  the 
location  of  the  power  deposition  and  the  activation  energy,  to  identify  “...the 
interrelationship  between  energy  transfer  occurring  at  a  high  rate  in  a  compressible 
medium  and  the  concomitant  motion  set  up  in  the  this  medium''.  Although  this  initial 
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value  study  is  limited  to  one-dimension  it  provides  useful  insights  into  the  impact  of 
localized  chemical  power  deposition  levels  on  induced  gasdynamic  processes. 

The  late  time  solutions  described  in  SKH  are  reexamined  to  ascertain  the  source  of 
relatively  small  solution  oscillations  that  appear  as  the  overdriven  detonation,  relaxes  to 
the  appropriate  stable  C.J,  state.  Optimization  of  the  Courant  number  value  and  the  grid 
size  used  in  the  MacCormack  method  are  used  to  demonstrate  that  the  oscillation 
amplitude  can  be  reduced  by  decreasing  the  mesh  size.  These  results  imply  that  the 
oscillation  has  a  numerical  origin  and  can  be  made  as  small  as  desired  by  grid  size 
reduction.  For  sufficiently  small  grid  size,  the  oscillations  in  several  sensitive  variables 
are  nearly  negligible. 

Prior  to  the  computational  study,  the  describing  equations  are  used  as  the  basis  for  a 
qualitative  discussion  of  the  physicochemical  processes  occurring  in  the  system.  A  novel 
description  is  given  of  the  relationship  between  localized  thermal  power  deposition  in  the 
gas,  the  appearance  of  confined  high  pressure  and  temperature  regions,  the  impact  of 
shock  waves  evolving  in  a  reactive  atmosphere,  and  the  gas  velocity  field  induced  by  the 
energetics  (Kassoy,  Kuehn,  Nabity,  and  Clarke,  2005).  It  is  also  shown  that  a 
transformation  of  the  variables  can  remove  the  activation  energy  from  the  Euler 
equations.  The  canonical  equations  are  cited  in  Eulerian  and  Lagrangian  form.  Further 
analysis  will  appear  in  a  future  paper. 

Finally,  computational  results  are  used  to  quantify  the  impact  of  thermal  power 
deposition  location  and  amplitude  as  well  as  activation  energy  on  the  evolution  of 
detonation  formation.  In  effect,  one-dimensional  models  are  usually  based  on  a 
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symmetric  initiating  disturbance  relative  to  the  geometrical  configuration  of  the  spatial 
region.  The  current  work  features  spatially  distributed  heating  on  the  microsecond  time 
scale  into  a  localized  region  at  a  finite  distance  from  a  confining  planar  boundary  (wall). 
The  initial  burst  of  thermal  power  into  the  asymmetrically  located  target  volume  causes  a 
significant  local  temperature  and  pressure  rise.  As  a  result,  strong  compressive 
gasdynamic  waves  propagate  both  upstream  and  downstream.  The  upstream  waves 
reflect  from  the  wall  and  evolve  into  strong  shocks  propagating  through  previously 
compressed,  partially  burned  gas  mixture.  The  wave  reflections  alter  the  overall  process 
of  detonation  formation  relative  to  that  found  with  symmetric  power  deposition. 

2  Formulation 

The  nondimensional  mathematical  model  for  a  compressible  gas  with  a  one  step 
exothermic  Arrhenius  reaction,  derived  by  SKIT  for  a  wide  range  of  physical  and 


chemical  conditions  as  well  as  length  and  time  scales  can  be  written  in  the  form: 

p  =  pT  (1) 

P,  +  {P“\  =°  ,  (2) 

p{ux+uu:)  =  -pjy  ,  (3) 

{pCj(y - J))(7;  +uT:)  =  -pu:  +  pQ  +  (pBqY/(y - l))exp(- 1/ s‘T )  ,  (4) 

K,  +uK  =-(5K/y)exp(-l/e*7')  ,  (5) 


The  thermodynamic  variables  (p,  p,  T)  are  nondimensional ized  with  respect  to  an 
undisturbed  initial  state  (p'f,  p'n,  T’t)  and  the  speed  u  with  the  undisturbed  speed  of  sound 
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o',  whereof  =/R'T'>.  Also,  the  inverse  activation  energy  parameter  e*  is  defined  with 
respect  to  T‘r 

The  nondimensionai  thermal  power  source  per  unit  mass  Q  is  absorbed  on  a 
dimensional  length  scale  characterized  by  C  with  a  duration  characterized  by  the  acoustic 
time  of  the  heated  domain^  =  tla[  It  follows  that 

a-lr/ir-WMle/W]  .  i-rtBjd.  r  =  C'JCl  (6) 

where  B'  is  the  preexponential  factor  in  the  Arrenhius  rate  law.  Further,  the  space  and 
time  variables  are  nondimension lized  by  using  C  and  C'la'n ,  respectively.  The 
nondimensionai  heat  of  reaction,  q  =  q'IC'pJ{ j',  where  q'  is  the  dimensional  value.  The 

remaining  variables  are  defined  in  standard  ways. 

Relative  to  the  terms  shown  in  Equations  (3)-(5),  SKI  I  also  show  that  transport  terms 
are  proportional  to  the  Knudsen  number  which  is  typically  very  small,  In 

general,  the  Euler-forms  of  the  equations  are  valid  when  ('  is  significantly  larger  than  the 
molecular  mean  free  path  and  far  smaller  than  the  axial  length  of  a  laboratory  detonation 
lube.  The  latter  condition  ensures  that  a  wave  process  initiated  near  the  left-hand  wall 
will  not  reach  a  right  hand  boundary  for  x-values  of  interest.  Further,  the  equations 
describe  the  time  history  of  a  spatially  varying  process  everywhere,  except  in  shocks 
which  are  described  by  jump  conditions  across  a  discontinuity. 

Initially,  the  undisturbed  system  is  described  by 

s  =  0  ,  p  =  p  =  T  =  Y  =  1,h  =  0  ,  2  >  0  .  (7) 


it 


Conditions  at  the  impermeable  left  hand  boundary  are  described  by 

2  =  0  ,  u  =  Y;  =  T,  =  0  .  (8) 

The  entire  reactive  event  is  initiated  at  s  -  0+  by  depositing  a  finite  amount  of  power  Q 
for  a  specific  period  of  time  into  a  localized  asymmetric  region  adjacent  to,  or  near  the 
left  hand  boundary,  defined  in  this  work  by: 

pQ  =  0  ,  0  <  z  <  z,  -  2  (9) 


=  3yf{s)  COS 


(zJ-2)<z<(z,+2) 


=  0 


z  >  z,  +  2 


(10) 

(11) 


where  z,  =2,  3  or  4  and 

/<.y)  =  0.7  [tanh  5  (5  —  £)-  tanh  5  (5-  1 0)j  (12) 

In  effect  the  initiating  power  deposition  terminates  very  close  to  the  switch  off  time 
A' =  10.  It  is  noted  that  the  spatial  structure  in  (9)-{l  i)  differs  from  that  employed  by  SKIT 
The  wave  system  evolving  from  the  initiating  power  deposition  process  propagates 
into  a  relatively  cold  gas  with  a  characteristic  spatially  homogeneous  thermal  explosion 

time  Ar  =e‘elu  /Bq,  which  is  typically  significantly  larger  than  the  time  scale  required 
for  a  CJ.  detonation  to  evolve  from  the  initiating  event  (e.g.,  a  =  0(10)).  For  example, 
the  parameters  used  in  SKI  I  (e*  =  0.08.  B=  15,  q  =  6)  give  se  =238,5.  It  follows  that  the 

change  in  T  in  front  of  the  oncoming  wave  system  will  be  very  small  (perhaps  1%) 
compared  to  that  in  the  system  for  x-values  of  interest.  Nonetheless,  it  is  useful  to  allow 
the  slow,  spatially  homogeneous  reaction  ahead  of  the  wave  system  to  occur  in  order  to 
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carry  out  the  smoothest  possible  initial  value  calculation.  In  effect,  the  lead  shock 
propagates  into  a  slightly  depleted  (perhaps  1%)  reactive  mixture. 


3  Properties  of  the  Describing  Equations 

A  qualitative  understanding  of  the  physical  processes  occurring  in  a  transient 
compressible  flow  with  localized  heat  addition  from  exothermic  reactions  can  be 
obtained  from  the  mathematical  properties  of  Equations  ( 1  )-(5).  In  particular,  it  is  useful 
to  demonstrate  the  relationship  between  thermal  power  deposition  into  the  gas  and  the 
resulting  induced  gasdynamic  processes,  including  the  appearance  and  evolution  of  shock 
waves. 

Equations  ( I ),  (2),  and  (4)  can  be  combined  [Kassoy  et  al.  2005]  to  show  that 


(y~0  jL+AL  e~'h'T  Cu  Dp 


(13) 


Y  CpT  /CpT  yCpp  Ds 


in  any  How  region  where  the  derivatives  are  defined  (e.g.,  not  including  shock 
discontinuities).  Here  D/Ds  =  (d/ds  +  ud/dz).  At  a  given  time  value  s,  one  may  integrate 
Equation  (13)  continuously  from  a  left  hand  boundary  (which  may  be  just  to  the  right  of  a 
shock  wave)  to  any  larger  value  z  prior  to  the  location  of  the  next  shock  downstream.  If, 
in  addition,  one  accounts  for  the  velocity  jumps  (Au)t  =u*  across  each  of  n  shock 
discontinuities  located  atz„; 
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where  the  integral  is  written  for  piecewise  continuous  functions  in  each  interval 
z„  <z<z1(J+)),  and  the  first  of  (8)  is  satisfied.  Here  {An),  may  be  positive  or  negative, 

depending  on  the  difference  in  velocity  just  to  the  right  of  (+)  and  just  to  the  left  of  (-),  a 
shock  which  may  move  to  the  right  or  left,  if  the  lead  shock  is  located  at  z,<„+1)(.s),  then 
the  velocity  jump  across  the  shock  must  be  equal  to  the  negative  of  the  induced  speed  just 
behind  the  shock:  (z1m)»+i  =—  w(zf,n+1))s)  s  — u~.  Here,  w"  is  the  “piston  speed"  supporting 

the  lead  shock  [Liepman  and  Roshko  1984b], 

Equation  (14)  can  be  used  as  a  diagnostic  tool  to  interpret  the  role  of  certain  physical 
effects  in  the  evolving  velocity  field.  The  first  term  in  the  integrand  represents  the  impact 
of  the  initiating  power  absorption.  In  general,  contributions  from  this  term  are  limited  to 
finite  space  and  time  domains.  In  the  specific  case  defined  in  (9)-(!2),  the  spatial  integral 
is  nonzero  for  (2,  -2)  <  2  <  (2,  +2)  and  s  <  1 0. 

The  second  term  in  the  integrand  represents  the  effect  of  chemical  heat  release 
relative  to  the  enthalpy  of  the  gas.  In  general.  0(1)  contributions  to  this  integral  are 
highly  localized,  existing  in  regions  of  significant  reactant  concentration  (Y  >  0),  and 
relatively  large  temperature  because  e  is  typically  small  (although  not  necessarily  in  the 
asymptotic  sense  of  taking  a  mathematical  limit).  The  amplitude  of  the  temperature 
dependent  function  w  =  e'  l!,:T/CpT  grows  relatively  rapidly  with  T  for  T  <()(\ie  )  and  the 
relative  increase  at  a  given  temperature  value  wjlw  ~  0(\le1^),  but  decreases  with 
increasing  temperature  levels.  Table  1  provides  values  for  w(7)  when  Cp  —  1  for  l  <  T< 
Me' .  In  hot  regions,  w  is  typically  10-2  or  larger  in  magnitude.  In  regions  where  the 
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temperature  is  relatively  low,  usually  prior  to  passage  of  a  sufficiently  strong  shock  wave, 
or  when  the  reactant  has  been  expended  Y  =  0  in  fully  burned  regions,  the  integral  of  the 
second  term  makes  little  or  no  contribution.  The  scaling  transformation  f-c'T  can  be 
used  in  the  second  term  to  show  that  then  Y  >  0,  the  integral  makes  a  sizable  contribution 
if  Bcjs'  =  0(1)  in  the  qualitative  sense. 


Table  1.  The  values  of  w(T)  for  several  activation  energy 
parameter  and  temperature  values  of  interest. 


iv  =  (l/r)e^*'r 

T 

£  =  .07 

£■*=.08 

II 

o 

1 

6.25  x  i  0"7 

3.73  x  I0-6 

4.54  x  10'5 

2 

3.95  x  I0-4 

9.75  x  I0‘4 

3.37  x  10"3 

3 

2.85  x  10‘3 

5.17  x  10"3 

1.19  x  I0-2 

5 

1.15  x  I0-3 

1.64  x  IQ“2 

2.71  x  10-2 

8 

2.10  x  10‘2 

2.62  x  10“2 

3.85  x  1 0-2 

Me 

2.58  x  10“2 

2.94  x  10“2 

3.68  x  10"2 

Fluid  particle  compression  is  described  by  the  last  term  in  the  integral.  In  local 
regions  of  partial  inertial  confinement  (nearly  constant  volume  combustion),  where  the 
relative  compression  ( 1  /p)Dp/Ds  is  important,  this  term  can  reduce  the  rightward  induced 
fluid  speed  or  reverse  it.  Thus,  some  thermal  power  deposition  results  in  an  induced  How 
field,  but  in  regions  of  near  inertial  confinement  a  pressure  rise  is  experienced,  and  little 
thermal  energy  conversion  to  kinetic  energy  (the  velocity  field)  takes  place.  Generally, 
inertial  confinement  occurs  in  regions  where  the  local  heat  addition  time  is  small 
compared  to  the  local  acoustic  time. 
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The  compression  term  can  also  be  interpreted  in  a  different  way  by  examining  the 
enthalpy  form  of  the  energy  equation  derived  from  (1),  (2),  and  (4) 

Cp  DT  (r-l)  l  Dp  f(r-l)\Q  t  BqY 
T  Ds  y  p  Ds  y  T  (y-\)T 

valid  in  regions  where  the  derivatives  exist.  The  first  term  on  the  right  hand  side 
provides  compression  heating  of  a  fluid  particle  as  a  supplement  to  the  other  thermal 
source  terms.  Equation  { 1 5)  also  provides  a  measure  of  the  deviation  of  the  fluid  particle 
processes  from  isentropic  conditions  that  would  prevail  if  the  thermal  source  terms  were 

absent.  In  addition,  one  can  observe  that  in  spatial  domains  where  £>  =  0,  where  the 
reactant  is  burned.  Y  —  0  and/or  where  the  temperature  is  relatively  low  (see  Table  ! )  the 
thermal  source  term  is  negligible.  Each  fluid  particle  behaves  isentropically,  although  the 
particle  entropy  will  vary  from  particle  to  particle,  because  the  physicochemical 
processes  occurring  in  any  particle  evolve  to  the  isentropic  condition  at  an  entropy  level 
that  depends  on  the  detailed  time  history  of  the  specific  particle.  Of  course,  even  in  the 
absence  of  source  terms,  the  entropy  of  a  fluid  particle  may  be  altered  suddenly,  due  to 
the  passage  of  a  shock  wave. 

Equation  (13)  also  implies  that  fluid  particle  compression  is  driven  by  thermal  power 
deposition  in  regions  of  small  or  zero  velocity  gradient.  In  such  a  region.  Equation  (2) 
can  be  used  to  show  that  Dp/Ds  *  0.  implying  that  a  nearly  constant  volume  process  is 
taking  place  during  a  brief  period  of  inertial  confinement. 

Significant  localized  increases  in  pressure  arc  one  type  of  gas  response  to  localized 
chemical  power  deposition.  The  appearance  of  a  high  pressure  domain  of  finite  extent  is 
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followed  by  a  relaxation  process  in  which  steep  fronted  compression  waves  propagate 
into  the  neighboring,  usually  moving  material.  Nonlinear  processes  cause  w-ave 
steepening  and  shock  formation.  Once  created,  shocks  may  be  located  instantaneously  at 
various  positions  z,„(.v),  mentioned  earlier. 

The  basic  describing  equations  in  Equations  (l)-(5)  can  be  put  into  Lagrangian  form 
by  using  the  transformation 

m(z,s)  =  Jf pdz  ,  (16) 


where  the  integral  is  taken  at  a  specified  time  value  s  in  a  piecewise  manner  to  account 
for  jumps  in  density  across  shock  discontinuities.  It  follows  that: 


C  T 

(r-. 0 


pu  =  T 

(17) 

25 

II 

(18) 

1 

II 

* 

5! 

(19) 

--PU.+Q+ 

(20) 

y. 

(21) 

where  u  =  1  ip  is  the  specific  volume.  The  analogue  to  (15),  for  the  case  of  constant 

specific  heats  CU=CP=  I ,  can  be  combined  with  (2 1 )  to  find 

-1)  (y-  1 

— ,-vpf-qYs  +  ^~ — 

7  7 

In  regions  unaffected  by  Q  (sufficiently  large  time  and/or  distance  from  a  thermal  power 
source),  one  can  integrate  this  expression  for  a  fixed  fluid  particle  (m  -  constant)  from  .v 
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(22) 
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=  0.  It  must  be  recalled  that  shocks  can  intersect  the  particle  at  a  sequence  of  times  s,.  At 
those  times  the  derivatives  in  (22)  are  undefined,  and  one  must  account  for  jumps  in  the 
thermodynamic  variables,  f  lowever,  no  discontinuity  in  Y  is  expected  at  s,  because  the 
chemical  response  to  a  change  in  state  occurs  on  a  finite  time  scale  following  shock 
passage.  It  follows  that 

r  =  I  +  g(l-r)+i^i  ["v^ds^AT,  ,  s„<s<s„,x  (23) 

f  Jo  & 

where  AT,  =  7^-7^“  and  the  +/-  superscripts  refer  to  locations  just  ahead  of  and  just 

behind  the  shock,  respectively.  Equation  (23)  shows  that  once  fuel  is  fully  consumed  Y  = 
0  in  s„  <  s  <  s*+],  the  temperature  is  a  superposition  of  the  adiabatic  flame  value  !  +  q ,  the 
accumulated  impact  of  the  compression  heating  time-history  described  by  the  integral 
expression  taken  in  piecewise  form,  and  the  net  effect  of  all  shocks  that  have  intersected 
the  particle  in  the  past.  It  is  noted  that  significant  net  compression  heating  and/or  the 
passage  of  strong  shocks  may  raise  the  particle  temperature  significantly  above  the 
adiabatic  value  1  +  q. 

The  Arrenhius  rate  law  embedded  in  the  energy  equation  in  Equation  (20)  can  be 
used  to  define  the  relative  local  chemical  reaction  time 

sCH  =\/Be~'lt'T  ,  (24) 

with  respect  to  the  duration  of  the  initiating  power  pulse,  s  —  0{ 1 ). 

Table  2  gives  values  o f  sen  for  a  small  temperature  range  of  interest  when  B  =  15 ,  the 
value  used  in  SKH.  As  long  as  reasonable  amounts  of  fuel  are  available,  q  =  (7(1)  and  Cv 
=  0(1),  then  Equation  (20)  can  be  interpreted  to  mean  the  substantial  changes  in  the  fluid 
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particle  temperature  will  occur  if  sch  is  not  too  large.  For  example,  the  numbers  in  Table 
2  suggest  that  during  the  early  phases  of  the  initiating  power  deposition  process,  when 
Q  =  G(  1}  and  T  =*  1+,  the  gasdynamic  phenomena  will  be  essentially  inert.  However, 
should  a  shock  wave  intersect  an  initially  cold  fluid  particle  and  raise  its  temperature  by 
as  little  as  a  factor  of  two,  the  value  of  sen  is  reduced  drastically  for  the  £  values  shown. 
It  follows  that  the  rate  of  heat  addition  from  chemical  processes  will  increase 
significantly  after  shock  passage  and  promote  significant  increases  in  fluid  particle 

temperatures.  The  accelerated  chemical  heating  process  is  itself  enhanced  by  increasing 

# 

£  , 


Table  2.  The  value  of  sch(T)  for  several  activation 
energy  parameter  and  temperature  values  of  interest. 


-  </**' 

T 

£  =  0.07 

£  =  0.08 

ii  ; 

p 

1 

1.07  x  10s 

1.79  x  I04 

1.47  x  103 

2 

84.5 

34.53 

9.89 

3 

7.8 

4.30 

1.87 

The  chemical  kinetics  term  in  (4)  can  be  used  to  define  a  nondimensional  local  (f) 


heat  addition  expression, 


(Lki 


{s,z)  = 


qB 


(rpyu/e')  Hr- 0 


pYe 


,-1/t  T 


(25) 
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used  to  ascertain  the  instantaneous  spatial  distribution  of  the  chemical  heat  release.  The 
global  instantaneous  heat  release  is  obtained  from  (25)  by  integrating  from  the  boundary 
to  the  lead  shock  location  z,  in  the  piecewise  sense. 


P 

?,»W  =  -r-A  \  ■  (26) 

nr- 1)  l 

Although  heat  is  released  in  front  of  the  shock  location,  the  magnitude  is  relatively 
insignificant  because  the  exponential  term  is  very  small,  as  in  Table  I.  Equations  (25) 
and  (26)  will  be  used  to  interpret  the  numerical  results  in  Section  6. 

4  Scaling  transformations 

The  relative  simplicity  of  (17)-(21)  has  enabled  the  discovery'  of  a  nonlinear 
transformation  set 


H  =  eVc'r  ,  p  =  Pfsm  .  v=V/ein  ,  u  =  U/e',n  (27) 

which  remove  the  dependence  on  £  from  (l)-(5)  when  qz  =0(1)  and  Qt :*  =0(1).  The 
results  arc 


PV  In  //  =  1 


(28) 


K=Um  (29) 

V,  =-Pm  (30) 


c„  H, 

(r-0  H(\nH)2 


-PV. 


If 


(31) 


Y  = 


BY 

YH 


(32) 
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Finally,  the  inverse  of  the  Lagrangian  transformation  in  (16)  can  be  written  as 

m 

z  =  s'~vl  jv (s,m)dm  (33) 

o 

The  transformation  in  (27)  desensitizes  the  describing  equations  (27)-(32)  with  respect  to 
the  value  of  s .  However,  it  is  recognized  that  the  activation  energy  parameter  value 
plays  a  role  via  boundary  conditions  of  the  lead  shock. 

Figure  (33)  also  implies  that  the  length  scale  transformation, 

z  =  X/e'12  (34) 

can  be  used  with  (25)  in  (l)-(5)  to  find  a  primitive  variable  canonical  equation  system; 

P  =  R/\nH  (35) 

R,+(RU)X=  0  (36) 

R{U,  +  UUx)  —  Pxfr  (37) 

-C„  ( H,  +  UHX )  =  - (r - 1 )  //  In  HUX  +  (r - 1 ) H  In2  H (*' 'Q )  +  In' 2  HB(qe  )Y  (38) 

Ys+UYx=-(by!Yh)  (39) 

where  R=  MV(p=emR). 

Finally,  the  global  heat  addition  function  in  (26)  can  be  transformed  to 

(?.,(j)/?)  =  («A(r-i))  I  (Rr/n)dx  ,  (40) 

0 

again  independent  of  e. 

Equation  (34)  implies  that  a  new  length  scale,  associated  w-ith  the  activation  energy 
parameter,  can  be  defined  by 
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x'  =  (rle'")x  , 


suggesting  that  if  the  initiating  power  deposition  occurs  on  the  C  -  length  scale,  then  the 
interesting  reactive  gasdynamics  occurs  on  the  scale  C’js  1  \  Numerical  solutions  for 
(28)-(32)  and  (35)-(39)  will  be  discussed  in  a  future  paper  to  ascertain  the  significance  of 
the  newly  discovered  transformation  in  (27). 

5  Numerical  Method 

The  numerical  solution  for  (l)-(5)  is  developed  by  using  a  second-order  accurate 
(space  and  time),  explicit  finite-difference  method  of  MacCormack  (1969),  including  a 
flux -corrected  transport  technique  described  by  Book  et  al  (1975).  Considerable  detail  is 
provided  in  SKH  and  in  Sileetn  (1989). 

The  CFL-stability  criterion  must  be  satisfied  at  all  mesh  points.  Operationally,  the 
time  step  must  be  found  from 

As  =  AA zj  ( |»|  +  a  )  ,  (41) 

where  A  is  the  Courant  number,  Az  is  the  selected  grid  size,  and  a  is  the  local  sound 
speed.  Extensive  testing  has  been  employed  to  find  optimized  values  of  A  and  A z,  as 
described  in  Section  6. 

As  an  example,  the  calculation  for  €  =  0.08,  B  =  15 ,  q  =  6,  y=  1 .4,  Ci/=  Cp  =  I  ,Pr  = 
0.72.  described  in  SKH  for  symmetric  power  addition,  has  been  redone  to  assess  the  role 
of  A  and  Az  on  the  time-history  of  the  solution.  Particular  attention  is  concentrated  on  the 
long-time  results  for  which  SKH  found  oscillatory  relaxation  of  an  initially  overdriven 
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detonation  wave  toward  an  average  CJ.  state  appropriate  to  the  system  parameters.  After 
considerable  mesh  refinement  at  K  =  0.98,  SKH  determined  that  at  Az  =  6.5  x  HP1,  the 
amplitudes  of  oscillations  relative  to  the  mean  values  where  typically  6.6%  for  the 
postshock  pressure.  8.1%  for  the  postshock  temperature,  and  2.5%  for  the  global  heat 
release  qch(s). 

Prior  to  the  appearance  of  oscillations,  the  transient  reactive  gasdynamics  appear  to 
be  relatively  insensitive  to  grid  spacing  for  6.5  x  10~J  <  Az  <  5  x  1 0-2.  This  is  not  too 
surprising,  given  that  the  most  rapid  transients  occur  on  a  timescale  of  about  five 
nondimensional  units  in  s  prior  to  detonation  formation.  Thus,  reductions  in  Az  in 
Liquation  (26)  cause  a  concomitant  decrease  in  the  time  step  which  is  already  small 
enough  at  Az  =  5  x  1CT2  to  resolve  the  time-history.  In  contrast,  the  oscillations  after 
detonation  formation  occur  on  a  time  scale  of  about  0.25  s  units,  so  that  a  reduction  of  Az 
by  a  factor  of  1 0  can  radically  improve  the  time  resolution. 

The  SKH  study  did  not  definitively  identify  physical  mechanisms  as  the  source  of 
the  oscillations.  Clearly,  numerical  processing  may  have  a  considerable  impact.  In 
particular,  the  large  absolute  temperature  and  pressure  values  behind  the  lead  shock  may 
promote  relatively  larger  errors  in  the  chemical  production  term  (see  (4))  just  behind  the 
shock,  arising  from  numerical  errors  in  the  temperature  itself.  This  is  particularly 
important  for  the  case  of  exponential  dependence  on  temperature.  For  example,  a 
specified  error  in  the  temperature,  6,  defined  by  T  =  Te  +  S,  where  Te  is  the  exact 


23 


solution,  creates  an  error  in  the  Arrhenius  term  F  ~  e  l/!fT  =  Fe  +  AF ,  where 
FF  =e  T,; ,  such  that 


dAF  F,S 


1 


-2 


(42) 


dr,  e't]  \_e'tk 

tt  follows  that  the  derivative  is  positive  (negative)  when  Te  <  M2e  for  a  prescribed 
constant  £  positive  (negative).  In  this  case,  the  absolute  value  of  the  error  AF  increases 
with  Te.  For  larger  Te,  the  error  decreases.  A  slightly  more  general  result  can  be  found  if 

S=  <%Tr). 


Improvements  in  the  solution  resolution  for  high  temperature  regions  in  which 
relatively  short  time  scale  processes  occur,  can  be  achieved  by  optimizing  the  values  of 
both  Az  and  A  relative  to  a  particularly  sensitive  solution  feature.  In  this  work,  the  time- 
history  of  the  pressure  behind  the  lead  shock  is  used  to  obtain  refined  values  of  the  two 
parameters,  and  the  results  are  supported  by  the  associated  behavior  of  qch(s)  defined  in 
(26). 

6  Optimization  Procedure 

Computational  solutions  for  0  <  5  <  40  have  been  obtained  for  the  thermal  power 
deposition  formula  in  SKH  (Equations  39-42),  with  £  =  0.08,  B  =  1 5 ,  q  =  6,  y=  1 .4,  Cv 
=  Cp  =  l,  for  a  variety  of  A  and  Az  values.  The  oscillation  amplitude  of  the  post  shock 
pressure  relative  to  the  mean  value  at  s  =  40  has  been  found.  First,  results  are  obtained 
for  Az  =  0.05  and  increasing  discrete  values  of  A,  beginning  with  A  =  0.8.  A  pronounced 
global  minimum  value  is  observed  at  A  =  0.91  as  seen  in  Fig.  I.  Thereafter,  results  are 
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obtained  for  a  sequence  of  increasingly  small  grid  sizes  down  to  Az  =  6.25  x  10~\  where 
the  relative  oscillation  amplitude  is  about  1%.  Subsequently,  results  are  obtained  at  two 
neighboring  values  of  X  =  0.91  to  show  that  away  from  the  minimum  very  slightly  larger 
oscillation  amplitudes  are  found. 

The  post-shock  pressure  as  a  function  of  time  s  is  given  in  Fig.  2a,b  for  a  relatively 
noisy  oscillatory  solution  (Az  -  0.05,  X  =  0,875)  and  a  more  accurate  result  (.42  =  6.25  x 
10_i,  X  =  0.92).  At  s  =  40  in  the  latter  the  mean  maximum  pressure  deviates  by  less  than 
1%  from  the  C.J.  value,  pcj-  35.73.  Larger  amplitude  oscillations  and  a  lower  average 
post-shock  pressure  are  observed  in  the  poorly  converged  case. 

In  Fig.  3a, b  results  are  given  for  the  time-history  of  the  global  heat  release  defined  in 
(26)  for  the  two  cases  in  the  previous  paragraph.  Here  again  there  is  a  significant 
reduction  in  the  amplitude  of  the  oscillation  for  the  optimized  case,  and  the  average  value 
is  within  a  fraction  of  a  percent  of  the  CJ.- value  of  83.1. 

These  results  suggest  quite  strongly  that  further  reductions  in  spatial  discretization  at 
an  optimized  value  of  X  would  produce  additional  reductions  in  the  oscillation 
amplitudes.  It  is  prudent  to  conclude  that  the  physical  long-time  relaxation  process  is 
smooth  rather  than  oscillatory,  and  that  the  small  oscillations  present  in  the  computational 
results  have  a  numerical  processing  origin.  In  this  regard  it  should  be  noted  that  crucial 
parameters  used  here,  q  =  6  and  e  ~  0.08,  differ  markedly  from  the  unstable  galloping 
detonation  wave  problem  studied  by  Fickett  and  Wood  (1966),  more  recently  by 
Bourlioux  et  al.  (1991)  and  many  others,  There  q  =  8.33,  e‘  -  0.02,  and  a  moving  piston 
supports  the  unstable  wave  propagation  process. 
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7  Computational  Results 


The  computational  results  are  to  be  given  a  thorough  physical  interpretation  in  order 
to  elucidate  the  fundamentally  transient  reactive  gasdynamies  occurring  as  a  planar  DDT 
evolves.  First,  a  description  is  given  of  events  following  power  deposition  into  a  gas 
volume  isolated  from  the  boundary,  where  z,  =  3  (refer  to  (10)).  Examination  of  the  time- 
history  of  the  spatial  distributions  of  the  dependent  variables  in  (l)-(5),  along  with  the 
local  and  global  heat  release  in  (25)  and  (26),  provides  insights  into  the  roles  of  shock 
wave  generation  and  multiple  reaction  center  (hot  spot)  formation  in  DDT  evolution. 
Subsequent  comparisons  with  results  for  z,  =  2  and  4  help  to  quantify  the  impact  of 
initiating  power  deposition  location  on  the  ensuing  reactive  gasdynamic  transients. 
Parametric  studies  for  the  case  z,  =  2  are  used  to  quantify  effects  of  the  initiating  power 
deposition  amplitude  and  the  activation  energy  of  the  one-step  reaction  on  the  evolution 
to  detonation. 

The  computational  results  can  be  interpreted  in  the  context  of  the  following 
dimensional  quantities; 

1 .  p\  =  I  atm,  T  =  300  K,  a,  =  350  mis 

2.  C  =  0.35  mm 

i.t'H  =  ('a  =\  f.is 

4.  q'  =  1.82x1  06J/kg,  5'  =  1.07x107/j,  y  =  1.4 

corresponding  to  the  nondimensional  parameters  B  -  1 5  and  q  =  6.  The  specific  heats  are 
assumed  to  be  constants  which  precludes  the  excitation  of  internal  modes.  In  addition, 
dissociation  is  not  included  in  the  model.  As  a  result,  the  temperature  extremes  are  larger 


26 


than  one  would  see  in  a  typical  reactive  gas  mixture  where  thermal  energy  can  be 
absorbed  by  excitation  of  rotation,  vibration  and  dissociation  of  molecules. 

Equations  (9-12)  can  be  used  to  find  the  dimensional  initialing  power  deposition  per 
unit  volume, 

£  * pQdz  =  1 .08 x  1 0"/ (**) w/m3  (43) 

where  f{s)  is  defined  in  (12).  The  integration  limits  in  (43)  imply  that  the  heated  region 
has  a  spatial  extent  of  1.4mm.  The  duration  of  the  deposition  is  10/x.t.  A  time  integration 
of  (43)  can  be  carried  out  to  find  the  total  energy  deposition  per  volume. 
E'd  =  1,53x10 V/m\  The  inert  temperature  rise  resulting  from  E'D  is  characterized  by 
0(1 0JK). 

The  numerical  results  have  been  obtained  with  a  uniform  spatial  resolution  of 
Ar  =  0.01 .  unless  otherwise  noted. 

7. 7  The  case  z,=  3,  e  =  0.0725 

The  temporal  variation  of  the  nondimensional  spatial  temperature  distribution  is 
given  in  Fig.  4  for  the  first  few  microseconds  of  the  deposition  process  when  z-t  =  3. 
Shortly  after  the  initiation  of  power  deposition,  lasting  to  s  =  10,  one  observes 
temperature  variation  beyond  the  heated  region,  1  <  z  <  5,  implying  that  the  an 
exothermic  chemical  reaction  has  been  initiated  in  the  neighboring  mixture,  the  result  of 
compression  wave  propagation  from  a  power  deposition  region.  Nondimensional  values 
of  the  temperature  greater  than  10  appear  for  s  =  2.5,  the  result  of  a  combination  of 
external  source  and  chemical  heat  release  (see  (23)).  In  fact,  species  distributions  in  Fig. 
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6  at  s  =  2  and  4  show  that  the  reactant  in  the  heated  region  is  totally  consumed  by  the 
latter  time  value. 

Thermal  power  absorption  and  exotherm icity  in  the  heated  region  (see  eq.  (14))  are 
the  sources  ol'gas  motion  and  compression  waves  that  propagate  both  to  the  left  and  the 
right.  The  properties  of  the  integrals  in  (14)  imply  the  induced  local  gas  Mach  number, 
defined  by  M ,  -  u'ja'  =  u]T 1,2 ,  is  less  than  one  for  the  source  deduced  in  (9)-(12)  when 
s  <  10.  Negative  velocity  values  (u)  are  possible  when  the  localized  fluid  particle 
1  DP 

compression  is  sufficiently  large.  Leftward  moving  gasdynamic  waves  reflect 

from  the  boundary  and  move  back  into  the  partially  or  fully  reacted,  heated  mixture. 
These  compression  waves  steepen  into  shocks  seen  in  the  lower  two  figures  in  Fig.  4.  The 
shock  Mach  number  in  the  latter,  based  on  data  at  s  -  3.4  and  3.5  is  Mt  w  1.3.  The 
observed  temperature  jump,  AT  *  3.4,  is  nearly  identical  to  the  expected  value  for  an  inert 
shock  when  y  =  1.4  (Liepman  and  Roshko  1957).  In  contrast  the  Mach  number  of  the 
evolving  lead  shock  moving  to  the  right  is  considerably  larger  because  the  sound  speed  of 
the  cold  gas  is  much  lower. 

A  reaction  zone,  represented  by  the  relatively  large,  but  smooth  temperature 
gradient,  trails  the  lead  shock,  it  moves  relatively  slowly  compared  to  the  rightward 
moving  shock  coming  from  the  left  seen  in  the  last  inset  of  Fig.  4  and  the  first  two  in  Fig. 
5.  A  coalescence  process  is  seen  in  the  latter,  leading  to  the  formation  of  a  more  sharply 
defined  reaction  zone  seen  in  the  third  and  fourth  insets  of  Fig.  5  for  s  ~  8.1  -  8.6.  The 
corresponding  species  result  can  be  seen  in  the  fourth  inset  of  Fig.  6  where  the  steep 
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concentration  gradients  arc  correlated  with  steep,  but  smooth  temperature  gradients 
referred  to  previously. 

The  large  temperature  values  observed  near  z  =  5  in  the  3rd  and  4th  insets  of  Fig.  5 
and  the  second  inset  in  Fig.  6  result  primarily  from  the  magnitude  of  the  externally 
supplied  thermal  energy  noted  below  (43).  Only  a  small  fraction  of  the  total  added 
energy  is  actually  needed  to  initiate  the  reaction  since  the  species  profiles  in  Fig.  6  show 
that  the  reactant  in  the  heated  region  is  burned  out  by  s  -  4,  while  the  external  heat 
addition  proceeds  until  s  =  10.  I.ater,  initiation  with  reduced  levels  of  external  thermal 
energy  is  considered. 

The  larger  lime  temperature  distribution  results  in  Fig.  5  and  those  that  follow  in  the 
first  two  insets  of  Fig.  6  show  that  about  I  Opts  after  the  heating  is  initiated,  a  combination 
of  reaction  processes  and  inert  gasdynamics  have  caused  a  shock-generated  temperature 
disturbance  out  to  about  r  =  22.  The  last  two  insets  in  Fig.  6  show  that  there  is  an 
un burned  pocket  of  mixture  adjacent  to  the  boundary.  To  the  right,  a  high  temperature 
region  of  combustion  products  is  found  (see  the  second  inset),  delimited  by  a  well 
defined  reaction  zone  (a  deflagration)  moving  into  a  gas  that  has  been  compressed  and 
heated  by  a  relatively  weak  lead  shock  (AT *  2). 

The  region  between  the  shock  and  the  reaction  zone  is  not  entirely  inert,  The  bumps 
in  the  temperature  distributions  represent  the  early  formation  of  localized  hot  spots 
moving  with  the  gas.  The  ripples  in  the  species  concentrations  at  s  =  6  -  10  (Fig.  6.  last 
inset)  represent  small  amounts  of  reactant  consumption  in  the  developing  hot  spots.  It  is 
noted  that  at  s  =  10  there  is  more  than  one  ripple  (and  bump)  in  the  curves,  representing 
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the  formation  of  more  than  one  localized  hot  spot.  This  result  may  be  contrasted  with 
that  found  in  SKH  where  symmetric  heating  leads  to  a  single  isolated  reaction  center.  In 
the  present  case  of  asymmetric  heating,  gasdynamics  w-aves  reflected  from  the  boundary 
are  the  presumptive  sources  of  the  multiple  spots. 

Further  evolution  of  two  localized  hot  spots  can  be  observed  in  the  first  inset  in  Fig. 
7  where  substantial  variations  in  the  temperature  can  be  seen  between  the  reaction  w-ave 
to  the  left  and  the  inert  shock  to  the  right.  The  corresponding  species  profiles  are  shown 
in  the  third  inset  of  Fig.  7,  where  the  reactant  consumption  pattern  is  clearly  discemable. 

Subsequently,  the  hot  spots  spread  and  coalesce,  produce  new  compression  waves 
moving  to  the  right  and  strengthen  the  lead  shock  by  the  lime  s  =  22  is  seen  in  the  second 
and  fourth  insets  in  Fig.  7.  At  $  =  1  8  one  observes  largely  unburned  pockets  of  reactant 
adjacent  to  the  wall  and  nearz  =  20.  In  addition  there  is  a  partially  burned  region  nearz  = 
27.  Beyond  that  there  is  a  smooth  reaction  front  in  the  region  30  <  z  <  37  which  separates 
fully  reacted  mixture  from  shocked,  unreacted  gas.  The  lead  shock  is  observed  to 
precede  the  reaction  zone  by  several  spatial  units. 

A  few  microseconds  later,  at  s  =  20  the  unburned  pockets  remain,  although  the 
second  of  the  two  has  moved  with  the  gas  a  few  spatial  units  to  the  left.  The  partially 
burned  region  now  located  near  z  =  28  has  nearly  vanished.  Meanwhile,  the  smooth 
reaction  front  has  moved  to  the  right  approximately  8  spatial  units.  The  right  side  of  the 
reaction  zone,  where  substantial  deviation  from  Y  =  1  exists  is  now  at  about  the  same  z 
location  as  the  lead  shock.  This  suggests  that  the  shock  has  been  sufficiently 
strengthened  to  enable  a  relatively  fast  reaction  to  exist  just  behind  the  discontinuity. 
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Finally,  at  s  -  22  one  can  observe  a  significant  temperature  bump  behind  the  lead 
shock,  representing  the  appearance  of  yet  another  localized  hot  spot  in  the  field.  The 
analogous  species  profile  in  the  fourth  inset  of  Fig.  7  shows  a  relatively  steep  reaction 
zone  just  behind  the  shock,  followed  by  a  partially  reacted  region  near  z  =  50  associated 
with  a  lower  temperature  seen  to  the  left  side  of  the  new  hump.  At  this  point  the 
detonation  wave  is  nearly  formed,  defined  as  a  self-sustained  wave  set  consisting  of  a 
shock  with  a  coupled  reaction  zone. 

Pressure  profiles  for  the  evolving  reactive  gasdynamic  process  are  shown  in  Fig.  8. 
The  first  inset  shows  a  localized  pressure  spot  at  s  =  2.  Further  pressure  enhancement  by 
heating,  subsequent  relaxation  (gas  expansion)  and  compression  wave  generation  to  the 
left  and  right  lead  to  the  shock  dominated  structure  at  s  =  4.  The  lead  shock  moves  to  the 
right,  followed  by  the  faster  shock  associated  with  wave  reflection  from  the  boundary. 
Coalescence  is  seen  in  the  second  inset  to  occur  just  after  s  —  6.  There  follows  a  period 
of  shock  decay  because  the  reaction  processes  described  above  are  decoupled  from  the 
lead  wave.  In  contrast  one  notes  lead  shock  strengthening  in  the  third  inset,  as 
compression  waves  generated  by  the  evolution  of  multiple  localized  hot  spots  propagate 
through  the  gas  toward  the  lead  wave.  Subsequently,  the  localized  chemical  heat 
deposition  just  behind  the  lead  wave  causes  a  very  rapid  enhancement  of  the  lead  shock 
and  one  sees  at  s  =  22  the  formation  of  the  pressure  profile  reminiscent  of  a  transient 
ZND  wave. 

The  time  history  of  the  local  chemical  heat  release  rate  distributions  w  ith  z  in  (25)  is 
given  in  Fig.  9  for  increments  of  time  between  19.6  and  26.  The  set  of  small  amplitude 
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spikes  to  the  left,  representing  heat  from  the  combustion  of  the  aforementioned  unbumed 
mixture  pockets,  are  considered  in  more  detail  in  Fig.  10.  The  largest  amplitude  profiles 
represent  the  consequences  of  rapid  chemical  reactions  just  behind  the  lead  shock.  The 
reaction  center  in  the  vicinity  of  z  -  40  at  s  =  19,6  evolves  into  a  bifurcated  structure  at 
s  =  20.6,  when  a  new  reaction  zone  appears  just  behind  the  strengthened  lead  shock.  The 
enhanced  maximum  values  reflect  the  increasing  post  shock  temperatures  as  the  lead 
shock  strengthens.  Further  evolution  of  this  structure  is  observed  at  s  =  21.6  and  23.0. 
The  shock-coupled  reaction  zone  at  the  latter  time  value  has  separated  from  the  nearly 
burned-out  reaction  center.  Subsequent  profiles  are  associated  with  an  over-driven, 
transient  detonation. 

Figure  10  shows  the  details  of  the  local  heat  release  distributions  for  the  unbumed 
pockets  at  four  time  values  (right  to  left)  beginning  with  s  -  19.6.  The  twin  local  peaks  at 
s  =  19.6  suggest  that  most  of  the  reaction  is  occurring  on  the  edges  of  the  pocket  Results 
at  subsequent  time  values  show  that  the  pockets  are  being  converted  to  the  left  by  the 
local  gas  flow.  The  profile  at  s  =  26  indicates  that  the  entire  pocket  has  begun  to  react, 
with  burnout  soon  to  occur. 

Species  profile  structural  features  in  Fig.  I  1  correlate  with  those  in  Fig.  9  and  help  to 
confirm  the  physical  interpretation  based  on  local  heat  release.  The  downward  spikes  in 
the  curves  at  s  -  20.6  and  21 .6  represent  the  birth  and  evolution  of  a  new  reaction  zone 
that  burns  out  at  larger  lime  values. 

Equation  (14)  is  used  to  produce  the  global  heat  release  rate  curve  in  Fig.  12  where 
data  are  recorded  at  intervals  of  As  =  0.0],  The  first  local  peak  is  associated  with  the 
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reaction  initiation  process  generated  by  thermal  power  deposition,  represented  by  the 
mesa-shaped  curve  beneath  the  peak.  A  subsequent  decline  in  chemical  power  output  is 
characteristic  of  the  early  phase  of  these  DDTs  where  close  coupling  of  shocks  and 
reaction  zones  is  absent.  In  the  vicinity  of  s  =  17,  one  observes  a  second  local  peak 
arising  from  the  appearance  of  the  multiple  reaction  centers  seen  in  Fig.  7.  The 
subsequent  global  maximum  is  correlated  with  the  appearance  of  the  multiple  reaction 
zones  at  s  =  20.6  and  21.6  in  Figs.  9  and  1 1.  The  rapid  decline  that  follows  occurs  when 
only  a  single  reaction  zone  is  coupled  to  the  shock  at  s  =  22.6.  the  result  of  a  sufficiently 
large  shock-induced  temperature  increase  which  facilitates  very  short  induction  times. 
The  fourth  local  maximum  arises  from  the  delayed  combustion  of  the  unburned  pockets 
portrayed  in  Fig.  10.  Finally,  the  global  heat  release  curve  asymptotes  (with  small 
oscillations  of  computational  origin)  to  the  steady  C.J.  value  of  qch  found  from  combining 
(14)  with  (5)  for  a  steady  state  wave, 
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The  corresponding  post-shock  pressure  value,  pej  —  35.64. 

Global  heat  release  rates  per  unit  area  in  a  shock  coupled  reaction  zone  (RZ)  can  be 
found  from  the  product  of  the  mass  input  rate  per  unit  area  and  the  heat  of  reaction  of  the 
reactant.  In  nondimensional  terms  (see  1-5), 
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<1bz=<1M,/{y-\)  , 

where  Ms  is  the  shock  Mach  number.  Eq.  (44)  is  found  for  a  CJ  wave.  Momentarily 
large  values  of  cjrz,  represented  by  the  spike  in  Fig.  12  imply  that  either  Ms  >  Mcj  or  that 
a  multiple  reaction  zone  structure  exists,  like  that  seen  in  Figs.  9  and  11.  In  fact,  the 
maximum  post-shock  pressure,  seen  in  Fig.  13 ,p  ~  40  is  compatible  with  a  shock  Mach 
number  Ms=  5.87,  which  corresponds  to  qxz  =  88.01,  only  about  half  of  the  value  found 
in  Fig.  12.  This  suggests  that  the  multiple  reaction  zone  structure  is  responsible  for  the 
transient  heat  release  spike. 

Equation  (25)  can  be  differentiated  with  respect  to  time  to  find  the  conditions  for  a 
global  maximum  in  qch  at  time  Sm.  It  follows  that, 

A^,(jM)exp[-lA*7;]  =  -  j  £-[/?Kexp(-l/fT)](/i  (46) 

where  fVs  is  the  shock  Mach  number  and  (ps>  Ts)  are  the  postshock  density  and 
temperature,  both  functions  of  W5.  The  integrand  represents  the  time  derivative  of  the 
rate  of  reactant  consumption.  It  is  noted  that  the  integral  must  be  negative  since  the 
lefthandside  is  positive.  It  follows  that  when  qck  is  maximized,  the  rate  of  chemical 
power  production  is  already  declining  in  the  average. 

The  time  dependence  of  the  local  pressure  maximum  is  given  in  Fig.  13  constructed 
by  identifying  the  largest  pressure  value  in  the  spatial  field  at  each  instant.  Lead  shock 
strengthening  generates  the  rapid  increase  in  pm ax  beginning  at  about  s  =  15.  The  global 
maximum  at  s  -  23.60,  associated  with  the  strongest  part  of  the  overdriven  detonation 
formation,  has  a  value  approximately  12.12%  greater  than  the  subsequent  CJ  value  of 
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35.64.  The  post-maximum  1.5%  “noise”  occurs  when  the  spatial  resolution  is  Az  =  0.0 1 . 
Further  reduction  in  grid  size  will  diminish  the  oscillations  as  the  results  in  Fig.  2  imply. 
The  mean  value  appears  to  be  asymptoting  to  the  C.J.  value  on  a  sufficiently  long  time 
scale.  One  should  note  that  the  maximum  power  level  qma*  occurs  at  s  =  21 .6.  before  that 
for  global  pressure  maximum  at  s  =  23.60.  This  sequence  of  results  is  characteristic  at 
many  different  parameter  values,  to  be  discussed  in  more  detail  in  later  sections.  This 
time  sequence  suggests  that  compression  waves  generated  by  the  localized  reaction  zone 
during  the  runup  to  the  heat  release  maximum  propagate  forward  to  the  lead  shock, 
maximizing  the  shock  pressure  a  few  acoustic  times  later. 

The  spatial  distributions  of  the  pressure  profiles  at  the  time  values  in  Figs.  9-1 1  arc 
shown  in  Fig,  14.  The  small  amplitude  disturbances  on  the  left  represent  relatively  weak 
compression  waves  propagating  toward  and  reflecting  from  the  planar  boundary.  The 
larger  profiles  to  the  right  show  the  shock  wave  strengthening  process  as  the  shock- 
coupled  reaction  zone  evolves  a  maximum  pressure  at  s  =  23.60  and  a  subsequent  decline 
as  depicted  in  Fig.  13. 

7.2  The  case  z,  -  4,  £  -  0.0725 

The  impact  of  moving  the  thermal  power  deposition  region  farther  away  from  the 
bounding  wall,  to  2  <  z  <  6,  can  be  observed  in  Figs.  15  and  16.  First,  it  is  noted  that  the 
general  features  of  the  global  heat  release  curve  are  quite  similar  to  those  in  Fig.  12  and 
that  the  initiation  power  curves  are  identical.  With  the  exception  of  the  initial  peak  in  qch 
at  about  s  =  3.  all  other  local  maxima  in  Fig.  15  occur  at  a  slightly  later  time  In  addition, 
each  local  maxima  is  larger  than  the  analogous  value  in  Fig.  12.  Finally,  the  shapes  of 


35 


the  absolute  maxima  curves  arc  nearly  identical,  although  they  occur  at  different  times 
and  have  different  absolute  values,  This  comparison  suggests  that  the  final  formation  of 
the  overdriven  deformation  is  independent  of  the  earlier  initiation  process. 

Subsequent  to  the  fourth  local  maximum  of  qCh  at  s  *  29.5,  due  to  the  heat  release 
from  the  consumption  of  the  previously  unbumed  pockets,  the  mean  value  of  the  curve 
asymptotes  to  the  C.J.  value  of  83. 1. 

A  comparison  of  results  in  Kigs.  13  and  16  indicates  that  the  absolute  maximum 
pressure  is  larger  and  occurs  later  when  Z/=  4.  Here  again,  it  is  observed  that  the  shapes 
of  the  primary  pmM  increase  with  time  are  nearly  identical  for  much  of  the  rise  to  the 
peak.  Finally,  the  pressure  peak,  associated  with  the  strongest  overdriven  detonation 
occurs  at  s  =  26,  later  than  that  for  the  absolute  gma*  at  23.8.  The  noisy  pressure 
relaxation  of  the  lead  shock  to  the  CJ.  value  appears  to  be  slower  than  that  of  the  global 
heat  release  in  Fig.  15.  Oscillatory  behavior  observed  in  Figs.  13  and  16  is  a  numerical 
artifact  related  to  issues  discussed  in  Section  4  and  described  in  Figs.  1  and  2. 

A  comparison  of  early  time  temperature  distributions  for  the  cases  z,  =  3  and  4  can  be 
found  in  Fig.  1 7a, b.  The  shocks  observed  in  both  cases  can  be  attributed  to  compression 
wave  reflection  from  the  boundary  surface  subsequent  to  the  explosion  of  reactant  heated 
by  the  initiating  power  deposition.  Shock  locations  for  z,  =  3  are  beyond  z  =  6  for  s  =  4.0, 
4.2  and  4.4.  In  contrast,  when  the  power  deposition  is  into  a  region  farther  from  the 
bounding  surface  (z,  =  4),  the  shocks  are  located  approximately  between  z  =  3  and  5.5.  It 
is  noted  that  the  shock  temperature  rises  observed  in  the  latter  case  are  within  a  few 
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percent  of  what  is  expected  for  a  y  =  i.4  shock  moving  at  the  observed  local  Mach 
number. 

7.3  The  case  z,  =  2. 

The  detonation  initiation  process  described  in  this  section  follows  thermal  power 
deposition,  defined  in  (9)-(12),  into  the  region  0  <  z  <  4.  This  placement  of  initiating 
power  minimizes  the  impact  of  shocks  reflected  from  the  bounding  surface,  in 
comparison  to  that  in  z,  =  3  and  4.  Results  are  presented  first  for  £  =  0.0725  and  varying 
power  deposition  amplitudes  with  the  goal  of  seeking  a  critical  value  below  which  the 
DDT  process  is  fundamentally  altered.  Secondly,  results  at  full  power  are  given  for 
several  values  of  s  to  determine  the  impact  of  activation  energy  on  the  detonation 
evolution  process. 

The  global  heat  release  variations  with  time  are  given  in  Figs.  I8{a-d)  for  power 
deposition  levels  of  100%,  80%,  60%  and  50%  when  e  =  0.0725.  The  full  power  result 
can  be  compared  with  that  in  Fig.  12  for  z,  =  3  to  ascertain  the  impact  of  moving  the 
center  of  power  deposition  closer  to  the  bounding  surface.  In  particular,  the  characteristic 
features  in  both  figures  are  quite  similar,  but  differ  in  timing  and  amplitude.  The  local 
maxima  occur  earlier  and  have  smaller  amplitudes  when  the  heated  region  is  closer  to  the 
bounding  surface,  A  summary  of  the  full  power  results  in  Figs.  1 2,  1 5  and  18a  is  given  in 
Fig.  19,  which  shows  how  the  amplitude  of  r/max  and  the  lime  to  qmax  varies  with  the 
location  of  the  center  point  of  the  heated  region  z,.  Both  appear  to  be  nearly  linear  for  the 
three  z, -values  considered.  Altered  gasdynamic  sequences  during  the  early  phases  of  the 
event  (see  Figs,  4  and  5  for  examples)  are  responsible  for  the  observed  differences. 
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The  sequence  in  Fig.  1  8  shows  that  decreases  in  the  power  deposition  amplitude  to  as 
small  as  50%  of  the  maximum,  delay  the  appearance  of  the  characteristic  features  and 
first  increases  and  then  decreases  their  amplitudes.  The  shape  of  the  £max  peak  is 
fundamentally  different  at  50%  of  full  power.  A  summary  of  q ma!(  amplitudes  and  times 
to  qm^  as  functions  of  initiating  power  deposition  is  given  in  Fig.  20  including  data  for 
many  input  values.  When  the  input  power  is  less  than  about  63%  of  maximum,  one 
observes  a  sharp  decline  in  qma*.  Concomitantly,  the  time  to  r/max  increases  smoothly  by 
more  than  100%  over  the  thermal  input  power  range. 

These  results  can  be  explained  by  recognizing  that  a  reduced  power  deposition  level 
leads  to  a  reduced  inert  temperature  rise  in  the  deposition  region  and  a  reduced  rate  of 
chemical  reaction  there.  The  maximum  value  of  the  initial  peak  in  the  global  heat  release 
declines  in  (a)-(d)  as  a  result.  The  increasingly  weak  gasdynamic  response  leads  to  a 
delay  in  the  detonation  formation  process. 

The  global  maximum  pressure  variation  with  time  is  given  in  Fig.  2!a-d  for  the  range 
of  power  input  levels  in  Fig.  I  8.  In  general,  the  absolute  maximum  pressure  (of  the  lead 
shock)  increases  as  the  power  input  level  decreases  while  the  time  to  that  value  increases. 
These  results  are  summarized  in  Fig.  22,  along  with  the  time  to  qmax.  There  appears  to  be 
a  change  in  character  for  the  pmax  time  when  thermal  power  input  is  less  than  about  55% 
of  the  maximum  value.  It  should  be  noted  that  the  time  to  qmax  precedes  the  lime  to  the 
absolute  maximum  pressure,  pm a*<abs)*  until  the  power  level  declines  to  about  57%  where 
a  crossover  occurs. 
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Local  heat  release  profiles  for  100%  power  deposition  are  given  in  Fig.  23  for  a 
sequence  of  limes  during  which  9max  occurs  at  s  -  18,8,  the  absolute  pressure  maximum, 
Pmaxubst  occurs  at  s  =  21.4  and  the  overdriven  detonation  is  born.  The  sequence  is  very 
similar  to  that  forz,  =  3  described  in  Fig.  9,  including  the  appearance  of  multiple  zones  of 
energy  release.  The  qm!a  occurs  when  the  original  reaction  center  and  the  shock-attached 
reaction  zone  are  both  generating  chemical  heat.  In  contrast,  /W{abs)  occurs  later  when 
the  only  substantive  heat  release  is  occurring  in  the  shock-coupled  reaction  zone 
associated  with  a  transient  overdriven  detonation. 

Species  concentration  profiles  in  Fig.  24  help  to  confirm  the  reaction  zone  structures 
seen  in  Fig.  23.  The  sequence  of  curves  on  the  right  side  show  how  a  new  reaction  zone 
appears  spontaneously  at  s  =  17.8.  The  trailing  reaction  zone  burns  out  at  s  *  22.4.  In 
Fig.  25,  the  pressure  distributions  with  z  are  given  for  the  same  time  sequence.  The 
emergence  of  an  unsteady  ZND  structure  is  observed,  driven  by  the  heat  release  from  the 
shock-coupled  reaction  zone. 

The  impact  of  activation  energy  on  the  detonation  initiation  process  can  be  observed 
in  Fig.  26a-d  where  the  global  heat  release.  </,  versus  time  is  given  for  £  =  0.0725, 
0.0750.  0.0775  and  0.0800  when  z,  =  2  and  the  initiation  power  is  100%.  A  decrease  in 
the  activation  energy  (increase  in  £  )  is  correlated  with  a  reduction  in  the  global  qmax  and 
its  occurrence  at  a  relatively  earlier  time. 

The  final  local  maximum,  associated  with  chemical  heat  release  from  the  combustion 
of  the  unburned  pocket,  is  observed  to  occur  closer  to  the  time  for  as  £  increases.  A 
significant  change  in  behavior  occurs  when  £  =  0.0775.  The  unbumed  pocket  is 
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consumed  during  ihe  period  of  most  intense  heat  release  and  detonation  formation.  As  a 
result,  the  fourth  local  maximum  disappears. 

The  value  of/;m:u((absj  in  Fig.  27a-d  declines  and  occurs  sooner  as  s  increases.  When 
s  <  0.08  the  absolute  maximum  pressure  exceeds  the  C.J.  value  and  pmax  approaches  the 
C.J.-value  from  above.  The  solution  properties  change  at  £  =  0.08  where  the  maximum, 
mean  shock  front  pressure  rises  monoton  ically  toward  the  C.J.  value  from  below.  The 
reduction  in  the  global  heat  release  maximum  seen  in  Fig.  26  is  compatible  with  reduced 
compression  wave  generation  and  lesser  front  shock  wave  enhancement. 

The  increase  in  qmax  with  decreasing  £*,  seen  in  Fig.  28  can  be  rationalized 
qualitatively  by  employing  the  thermal  gradient  concept  (spatial  gradient  in  the  induction 
time)  (Zeldovich  et.  al,  1970  and  Zeldovich,  1980).  A  reactive  mixture  fluid  particle 
compressed  and  heated  by  passage  of  the  lead  shock  releases  most  of  its  chemical  energy 
in  a  brief  burst  of  power  within  a  thin  reaction  zone  at  the  end  of  the  much  longer 
induction  zone,  the  result  of  the  exponential  dependence  of  reaction  rate  on  temperature 
(Kassoy,  1977).  Higher  activation  energy  mixtures  are  characterized  by  relatively  long 
induction  periods,  Extended  induction  periods  imply  that  a  relatively  large  mass  of 
nearly  unreacted  mixture  can  be  present  behind  the  shock  in  a  higher  activation  energy 
system,  During  a  period  of  shock  strengthening,  seen  in  Figs.  1 3  and  14,  the  induction 
time  for  a  sequence  of  neighboring  fluid  particles  is  actually  reduced,  leading  to 
possibility  of  multiple  reaction  zones  seen  in  Fig.  1 1.  Figure  28  also  shows  increasing 
values  of  post-shock  pressure,  /Jmax(abs)  with  decreasing  e*.  the  result  of  enhanced  local 
chemical  power  deposition  with  increases  in  activation  energy.  Concomitant  compression 
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wave  generation  leads  to  lead  shock  strength  enhancement.  Irregular  location  of  some 
numerical  data  points  arise  from  sensitivity  to  the  sampling  process. 

A  nearly  monotonic  linear  increase  of  the  time  to  7maX  with  decrease  in  £*  is 
observed  in  Fig.  29.  a  logical  outcome  of  the  lengthened  induction  time  with  increasing 
activation  energy.  The  time  to  /Wobo  has  a  minimum  near  e*  =  0.078.  with  a  very  large 
increase  for  further  decreases  in  s*.  This  can  be  interpreted  to  mean  that  when  the 
chemical  power  deposition  level  falls  below  a  critical  value  (see  Fig.  26  c  and  d).  the 
compression  wave  generation  process  cannot  enhance  the  lead  shock  strength  beyond  the 
C.J.  value  (see  Fig.  27d). 

8  Conclusions 

The  essential  properties  of  detonation  initiation  transients  are  described  for  stable, 
planar  detonation.  The  results  described  in  the  present  work  and  that  in  SKH  suggest  that 
the  transient  dynamics  of  stable,  planar  DDT  evolution  are: 

1.  dependent  on  the  characteristics  of  the  initial  power  absorption  process,  including 
the  location  relative  to  boundaries  present, 

2.  quite  sensitive  to  the  activation  energy  for  a  given  level  of  power  deposition, 

3.  characterized  by  the  spontaneous  appearance  of  several  localized  hot  spots 
(reaction  centers)  resulting  from  the  conditioning  of  isolated  volumes  of  reactive 
mixture  by  the  complex  compression  wave  structure  inherent  in  a  compressible 
system  with  thermal  power  addition  from  external  and  chemical  sources,  and  is 

4.  likely  to  lead  to  an  overdriven  detonation  which  relaxes  smoothly  to  a  C.J.  wave, 
unless  the  activation  energy,  E\  is  sufficiently  small.  The  overdrive  is  relatively 
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small  in  terms  of  the  maxi  mum  post-shock  pressure  and  lead  shock  Mach  number, 
although  the  global  heat  release  rate  peak  is  significantly  larger  than  the  C.J.  value 
for  a  steady  wave.  Multiple  reaction  zone  combustion,  facilitated  by  continuous 
reduction  in  the  induction  time  of  a  sequence  of  reactive  fluid  particles  passing 
through  a  continually  strengthening  shock  wave,  can  account  for  the  observed 
heat  release  dynamics.  For  example,  the  distance  traveled  from  the  shock  by  an 
«th  fluid  particle  during  the  induction  time  is  characterized  by  Z„  =  o-,„an  where 
crm  is  the  induction  time  for  that  particle  and  o„  is  a  measure  of  the  post  shock 
speed  relative  to  the  shock.  Both  quantities  are  declining  with  increasing  shock 
strength.  It  follows  that  the  location  of  large  heat  release  for  a  large  n  particle  can 
be  closer  to  the  shock  than  that  for  a  small  n  particle,  and  both  can  occur  at  the 
same  time. 

The  optimization  study  shows  that  the  small  oscillations  seen  in  the  SKH  results 
during  the  overdriven  detonation  relaxation  process  to  a  steady  C.J.  wave  can  be 
systematically  reduced  by  the  proper  choice  of  the  (Az,X)  combination.  The  optimized 
values  have  been  employed  for  the  asymmetric  power  deposition  DDT  results  described 
in  the  present  work. 
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Figures 


Deviation  of  Pressure  Prom  Mean  Maximum  Pressure  (e  -  0  08) 


l;igure  I .  Relative  post-shock  pressure  oscillation  amplitude  (%)  vs. 
Courant  number  (X)  with  respect  to  the  grid  size  (Az)  for  SKI  I 
parameters. 
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Figure  2.  The  post-shock  pressure  variation  with  time  for  a 
poorly  resolved  case  (a),  and  a  more  accurate  result  (b),  using 
SKH  data. 
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Figure  3.  The  global  heat  release  variation  with  time  for  a  poorly 
resolved  case  (a),  and  a  more  accurate  result  (b),  using  SKH  data. 


Figure  4.  The  spatial  temperature  distribution  For  several 
time  values  to  s  =  3.5  for  Zj  =  3  and  e*  =  0.0725. 


Figure  5.  The  spatial  temperature  distribution  for  several  time 
values  from  s=3.9  to  s  =  8.6  for  z,  =  3  and  £*  =  0.0725. 
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Figure  6.  The  spatial  temperature  and  species  distributions  for 
several  time  values  to  s  =  10  for  z,  =  3  and  e*  =  0.0725, 


Figure  7.  The  spatial  temperature  and  species  distributions  for 
several  time  values  from  s  =  12  to  s  =  22  for  z,  =  3  and  e*  = 
0.0725 


Figure  8.  The  spatial  pressure  distributions  for  several  time 
values  to  s  =  22  for  Zj  =  3  and  e*  =  0.0725. 
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igurc  9.  The  lime  history  of  the  local  chemical  heat  release 
slributions  with  z,  defined  in  Eq,  25,  for  the  time  values  s  =  19.6, 
).6,  21.6,  23,  24,  25  and  26  when  zp3  and  c*  =  0.0725. 


Figure  10.  The  time  history  of  the  local  chemical  heat  release 
distributions  with  z  in  the  unburned  mixture  pockets  for  the  time 
values  s  -  19.6,  21.6,  24  and  26  (right  to  left)  when  Zj  =  3  and  e*  = 

0.0725. 


Figure  1 1.  The  time  history  of  the  species  distribution  with  z  for  the 
time  values  s  =  19.6,  20.6,  21.6,  23,  24,  25,  and  26  when  Z\  =  3  and  t* 
=  0.0725.  The  seven  curves  beyond  z  =  25  are  read  left  to  right  with 
increasing  time.  The  unbumed  mixture  pockets  move  to  the  left  as 
time  increases. 
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Figure  12.  The  evolution  of  the  global  heat  release  with  time, 
defined  in  Eq.  26,  exhibiting  four  local  maxima  and  asymptoting  to 
the  CJ.  value  for  a  steady  detonation  wave  for  z,  =  3  and  e*  = 
0.0725.  The  mesa-like  curve  in  the  lower  left  hand  corner  represents 
the  initial  power  deposition.  The  absolute  maximum  global  heat 
release,  qmm  —  1 66,08,  occurs  at  s  =  21.6. 


Figure  13.  The  time  dependence  of  (he  instantaneous  global 
maximum  pressure,  p™,,,  corresponding  to  (he  appearance  of  an 
overdriven  detonation  and  subsequent  relaxation  to  the  CJ.  value  for 
a  steady  detonation  wave,  for  z,  =  3  and  e*  =  0.0725.  The  absolute 
maximum  pressure,  pmax  =  39.96,  occurs  at ■=  23.6. 


Figure  14.  The  time  history  of  the  pressure  distributions  with  z 
for  the  time  values  s  =  19.6,  20.6,  21.6,  23,  24,  25,  and  26  when 
Zj  =  3  and  e*  =  0.0725,  including  the  effect  of  relatively  weak 
compressions  reflecting  from  the  wall. 
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Figure  15.  The  evolution  of  the  global  heat  release  with  time, 
defined  in  Eq.  26,  when  z,=4  and  e*  =  0.0725.  The  absolute 
global  maximum  heat  addition,  gmax  =  176.61,  occurs  at  $=  23.8. 


Figure  16.  The  time  dependence  of  the  instantaneous  global 
maximum  pressure,  pmax,  corresponding  to  the  appearance  of  an 
overdriven  detonation  and  subsequent  relaxation  to  the  C.J.  value 
for  a  steady  detonation  wave  when  Z(=4  and  e*  =  0.0725,  The 
absolute  pressure  maximum,/?™,  -  40.78,  occurs  at  s *=  25.8. 


Figure  17.  The  time  history  of  the  spatial  temperature 
distributions  with  z  to  show  impact  of  the  initial  power 
deposition  location  when  e*  =  0.0725:  a.  Zj=3,  s=4.0,  4.2  and 
4.4;  b.  Zj=4,  s=4.0,  4.2  and  4.4. 
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Figure  18.  The  evolution  of  the  global  heat  release  with  time,  defined  in  Hq.  26, 
exhibiting  four  local  maxima  and  asymptoting  to  the  C.J.  value  for  a  steady 
detonation  wave  when  z,  =  2  and  e*  =  0.0725.  The  mesa-like  curve  in  the  lower 
left  hand  corner  represents  the  initial  power  deposition:  a.  100%  power 
deposition,  b.  80%  power  deposition,  c,  60%  power  deposition  and  d.  50% 
power  deposition.  The  corresponding  absolute  maximum  global  heat  addition 
values  are  qmm  =  158.87  at  s  =  18.8,  167.31  at  s  =  23.2,  181.66  at  s  —  33.6  and 
1 57.94  at  46.8,  respectively. 


Figure  19.  The  variation  of  the  absolute  maximum  global  heat 
release,  r/max,  and  the  time  to  qmm  with  z,  for  100%  power 
deposition  when  e*  =  0.0725. 
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Figure  20.  The  variation  of  the  absolute  maximum  global  heat 
release,  r/ma„  and  the  time  to  qm ax  as  a  function  of  the  %  of  initial 
power  deposition,  when  Zj  =  2  and  e*  =  0.0725. 


Figure  21.  The  time  dependence  of  the  instantaneous  global  maximum 
pressure,  corresponding  to  the  appearance  of  an  overdriven  detonation 
and  subsequent  relaxation  to  the  C.J.  value  for  a  steady  detonation  wave 
when  z,  =  2  and  e*  =  0.0725:  a.  100%  power  deposition,  b.  80%  power 
deposition,  c.  60%  power  deposition  and  d.  50%  power  deposition.  The 
corresponding  maximum  pressures  are  38.77  at  s  =  20.8,  40.74  at  s  =  24.8, 
46.91  at  s  =  34.6,  and  58.48  at  s  =  46.4. 
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Figure  22,  The  variation  of  the  absolute  pressure  maximum,  /Wtabsj, 
and  the  time  to  />ma*(abs)  for  a  range  of  initial  power  deposition  levels 
when  zr-  2  and  c*  =  0.0725. 
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Figure  23.  The  time  history  of  the  local  chemical  heat  release 
distributions  with  z,  defined  in  F,q,  25,  for  the  time  values  s  =  16.8, 
1  7.8,  1 8.8,  1 9.8,  21 .4,  22.4  and  23.4  when  Zj=  2  and  e*  =  0.0725. 


1 

09 

OS 

07 

OS 

OS' 

04 

0  3 

02 

01 

0 

0 

F 

t\ 

Zi 


r - t - t=t - : - 1 - 1 - 1 - 1 - r 


10  20  30  40  SO  60  70  SO  90  100 


z 


igure  24.  The  time  history  of  the  species  distributions  with  z  for 
e  time  values  s  -  16.8,  17.8,  18.8,  19.8,  21,4,  22.4  and  23.4  when 
=  2  and  e*  =  0.0725. 
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igure  25.  The  time  history  of  the  pressure  distributions  with  z  for 
e  time  values  s  =  16.8,  17.8,  18.8,  19.8,  21.4,  22.4  and  23.4  when 
=  2  and  e*  =  0.0725. 


Figure  26,  The  dependence  of  ihe  global  heat  release,  q,  with  lime,  s,  on  the 
activation  energy  parameter  E*:  a.  £*=0.0725,  £*=0.0750,  £*=0.0775  and 
£*=0.080  when  z,  =  2,  The  fourth  local  maximum  is  absorbed  into  the  right 
side  of  the  peak  associated  with  the  absolute  maximum  for  the  larger 
parameter  values.  The  corresponding  global  heat  addition  maxima  arc;  a.  qmax 
=  158.9  at  s  =  18.8,  b.  qm„  =  142.0  at  s  =  15.4,  c.  qmax  =  130.6  at  s  =  13.2,  and 
d.  qmax  =  1 24.0  at  s  =  11.4. 


Figure  27.  The  time  dependence  of  the  instantaneous  global  maximum 
pressure, pmax,  as  a  function  of  the  activation  energy  values  in  Fig.  26,  when  Zj- 
2  with  100%  power  deposition.  When  £*=0.080  the  pmax  variation  with  time 

differs  from  that  found  for  smaller  values  of  £  .  The  corresponding  absolute 
maximum  pressure  are  38.77  at  s  =  20.8,  36.99  at  s  =  1 7.4,  35.9668  at  s  =  1 5.2, 
and  35.29  at  5  =  32.0. 
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Figure  28.  The  variation  of  qm*x  (left  scale)  and  />max<abs>  (right 
scale)  with  activation  energy  parameter  values,  £*,  for  zt  -  2. 
Both  quantities  are  sensitive  to  small  variations  in  8** 
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Figure  29.  The  dependence  of  the  time  to  </max  and  /w<abs)  on  the 
activation  energy  parameter  value,  £*for  z,=  2.  The  former  shows 
a  nearly  linear  monotonic  decline  with  increasing  values  of  £*, 
while  the  latter  has  an  absolute  maximum  near  £*  -  0.0775. 
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